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Deformed shell structures in nuclear mean-field potentials are systematically investigated as functions of 
deformation and surface diffuseness. As the mean field model to investigate nuclear shell structures in wide 
range of mass numbers, we propose the radial power-low potential model, V (x r", which enables us a simple 
semiclassical analysis by the use of its scaling property. We find remarkable shell structures emerge at certain 
combinations of deformation and diffuseness parameters, and they are closely related to the periodic orbit bi- 
furcations. In particular, significant roles of the "bridge orbit bifurcations" for normal and superdeformed shell 
structures are pointed out. It is shown that the prolate-oblate asymmetry in deformed shell structures is clearly 
understood from the contribution of bridge orbit to the semiclassical level density. The roles of bridge orbit 
bifurcations to the emergence of superdeformed shell structures are also discussed. 

PACS numbers: 21.10-k, 03.65.Sq 



I. INTRODUCTION 

Shell structures in single-particle energy spectra play es- 
sential roles in nuclear ground-state deformations. As a phe- 
nomenological mean-field potential, modified oscillator (MO) 
and Woods-Saxon (WS) models are successfully employed in 
shell correction approach. For a simpler and qualitative de- 
scriptions of the properties of shell structures, harmonic oscil- 
lator (HO) and infinite-well (cavity) potentials are frequently 
aided for light and heavy systems, respectively. Axially- 
symmetric anisotropic HO potential model successfully ex- 
plain the magic numbers of light nuclei, emergence of su- 
perdeformed shell structures, etc. For heavier nuclei, radial 
profile of the potential around the nuclear surface becomes 
more sharp and it looks more similar to square well poten- 
tial. In order to prevent from complexity of treating contin- 
uum states, WS potential is sometimes approximated by in- 
finite well potential (cavity). The cavity system, as well as 
HO system, is integrable under spheroidal deformation due to 
the existence of a non-trivial dynamical symmetry, and sev- 
eral classical and quantum mechanical quantities are obtained 
analytically. It also accepts several useful techniques to calcu- 
late quantum eigenvalue spectra, since Schrodinger equation 
is equivalent to the Laplace equation with Dirichlet boundary 
condition. 

The HO and cavity systems have significant difference in 
deformed shell structures. The single particle spectra in de- 
formed HO potential are symmetric in prolate and oblate 
sides, while those in cavity system are apparently asymmet- 
ric. Such asymmetry has been considered as the origin of so 
called prolate-shape dominance in nuclear ground-state defor- 
mations; a well-known experimental fact that the most of the 
ground-states of medium-mass to heavy nuclei have prolate 
shapes rather than oblate shapes. This prolate-shape domi- 
nance has been reproduced theoretically in microscopic cal- 
culations. Tajima et al. made a systematic Nilsson-Strutinsky 
calculations throughout the nuclear chart jit] and examined 
the distribution of ground-state deformations. In order to 



pin down the essential parameter which causes prolate-shape 
dominance, they varied the strengths of l^ and Is terms in 
Nilsson Hamiltonian, and found that the prolate-shape domi- 
nance occurs under strong correlation between and Is terms. 
The recent analysis by Takahara et al. based on Woods- 
Saxon-Strutinsky calculations also support those resultsj^]. 
Hamamoto and Mottelson compared the oblate and prolate 
deformation energy in spheroidal HO and cavity models and 
found the prolate-shape dominance only in cavity model. 
They considered the origin of that prolate-shape dominance as 
the asymmetric way of level fanninngs in prolate and oblate 
sides for a potential with sharp surface, and show that the 
above asymmetry is explained from the different roles of in- 
teraction between single-particle levels in prolate and oblate 
sideslH. 

We expect that the semiclassical periodic-orbit theory 
(POT) holds the key for deeper understandings of above shell 
structures responsible for prolate-shape dominance. In semi- 
classical theory, quantum mechanical quantities can be related 
to classical mechanical variables associated with classical or- 
bits. In the periodic-orbit theory (POT), density of states is 
represented as the sum over contribution of classical periodic 
orbits: 

giE)=m+LME)cos^^-^^y (1.1) 

The sum on the right-hand side is taken over all classical pe- 
riodic orbits j3 which exist for energy E. Each periodic orbit 
(PO) j5 changes its shape with increasing energy E, and the 
action integral Sp = f^p ■ dr along the orbit is, in general, a 
monotonically increasing function of E. Thus, each cosine 
term in the PO sum dl.ll ) is a regularly oscillating function of 
energy whose period of oscillation 5E is given through the 
relation 

5S,^^^8E^2Kh, , (1.2) 

P dE 7|3 ' 

where = dSp /dE is time period of the orbit )3. Therefore, 
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a gross shell structure (large 5E) is associated with short pe- 
riodic orbits (small Tp). 

The relation between coarse-grained quantum level density 
oscillations and classical periodic orbits in spherical cavity 
model was first discussed by Balian and Bloch|4]. They show 
that the modulations in quantum level density oscillations are 
clearly understood as the interference effect of periodic orbits 
with different lengths. This idea has been successfully applied 
to the problem of supershell structure in metallic clustersoll. 
Strutinsky et al.llet] applied periodic orbit theory (POT)||lla] 
to the cavity model with spheroidal deformation and discussed 
the properties of deformed shell structures in medium-mass to 
heavy nuclei in terms of classical periodic orbits [|6t]. Frisk 
made more extensive POT calculations to reproduce quan- 
tum level density by the semiclassical formulaj'sl. He also 
suggested the relation between classical periodic orbits and 
prolate-oblate asymmetry in deformed shell structures, which 
might be responsible for prolate-shape dominance discussed 
above. Those works have proved a virtue of semiclassical 
POT for a clear understanding of properties of finite quantum 
systems. 

It should be emphasized here that unique deformed shell 
structures are developed when the contribution of certain peri- 
odic orbit are considerably enhanced. The magnitude of shell 
effect is related with the amplitude factor in Eq. (II. 11 1. This 
amplitude factor has important dependency on the stability of 
the orbit, which is generally very sensitive to the potential pa- 
rameters such as deformations. In particular, stability factor 
sometimes exhibit significant enhancement at periodic orbit 
bifurcation, where new periodic orbits emerge from an exist- 
ing periodic orbit. Near the bifurcation point, classical orbits 
surrounding the stable periodic orbit form a quasi-periodic 
family, which make coherent contribution to the level density. 
This is an important mechanism for the growth of deformed 
shell structures. 

A typical example is so-called superdeformed shell struc- 
ture. It is known that single-particle spectra exhibit remark- 
able shell effect at very large quadrupole-type deformation 
with axis ratio around 2:1. In anisotropic HO model, this shell 
structure is related with the periodic orbit condition; all the 
classical orbits become periodic at (0± = 2(0- and they make 
very large contribution to the level density fluctuation. In the 
cavity model, one also finds significant shell effect around 2:1 
deformation, and it is related with the bifurcations of equa- 
torial periodic orbits through which three-dimensional (3D) 
periodic orbits emerge Igl fioll . It should be interesting to ex- 
plore the intermediate situation between the above two limits, 
which might correspond to the actual nuclear situation. 

Our purpose is to understand the transition of deformed 
shell structure from light to heavy nuclei in terms of classi- 
cal periodic orbits. It requires a mean-field hke WS potential 
model. Periodic orbits in spherical and deformed WS poten- 
tials have been examined by Arvieu et al.[l 1, 12], but their re- 
lation to quantum level densities have not been discussed. The 
inner region of WS potential is nicely approximated by power- 
law potential which has simpler radial dependence V cx r". 
This approximation simplify both quantum and classical cal- 
culations and one can achieve clear quantum-classical corre- 



spondence via Fourier transformation technique. 

Thus, in the current paper (Part I), we focus on the radial de- 
pendence of the mean field potential (effect of surface diffuse- 
ness, described by the term in the Nilsson model) and ex- 
amine the shell structures systematically as functions of defor- 
mation and surface diffuseness. As pointed out by Tajima et 
al., spin-orbit coupling play also an important role in prolate- 
shape dominance. The effect of spin-orbit coupling will be 
discussed in the forthcoming paper (Part II). 

This paper is organized as follows. In Sect. |II] we dis- 
cuss the quantum and classical properties of power-law po- 
tential model. The scaling properties of the model are de- 
scribed and the Fourier transformation techniques are formu- 
lated. In Sect, [nil quantum mechanical density of states and 
shell structures in spherical power-law potential are examined. 
Some analytic expressions for periodic orbit bifurcations and 
semiclassical formula are given, and quantum-classical corre- 
spondence is discussed. It is shown that bifurcations of cir- 
cular orbit play essential role in unique shell structures which 
appear at several values of radial parameter a. In Sect. |IV] 
shell structures are examined against spheroidal deformation 
parameter Semiclassical origin of prolate-oblate asymmetry 
in deformed shell structures and that of superdeformed shell 
structure are investigated. Special attention is paid on what 
we call "bridge orbit bifurcations". Section |V] is devoted to 
summary and conclusion. 



II. THE POWER-LAW POTENTIAL MODEL 

A. Deflnition of the model 

It is known that the central part of the nuclear mean field 
potential is approximately given by the Woods-Saxon (WS) 
model 



VwsW = - 



w 



l+exp{{r-RA)/a}' 



(2.1) 



The depth of the potential is W ~ 50 MeV, surface diffuse- 
ness is fl ~ 0.5 fm and nuclear radius is ^ 1.2A'''^ fm 
for a nucleus with mass number A. For small A, the inner 
region of this potential can be approximated by harmonic os- 
cillator (HO). For large A, the inner region of the potential 
is flat {V « —W) and sharply drops to zero at the surface, 
which looks more similar to square well potential. In Ref. J^, 
shell energy of deformed WS potentials are compared with 
those for HO and infinite square-well (cavity) potentials. De- 
formed shell structures in WS model are similar to those of 
HO model for light nuclei, while they are more like those of 
cavity model for medium-mass to heavy nuclei. Our aim is 
to understand the above transition of deformed shell structure 
from the view point of quantum-classical correspondence. For 
the above purpose, we take the radial dependence of the po- 
tential as r", which smoothly connect HO (a = 2) and cavity 
(a = oo) potentials by varying the radial parameter a: 
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FIG. 1. Profiles of power-law potentials ( 12. 2t fitted to the Woods- 
Saxon potentials ( 12.1b for mass numbers A = 4, 60 and 240. Values 
of radial parameter a are determined by Eq. ( 12.4b . 



This power-law potential, having simple radial profile, is easy 
to treat in both quantum and classical mechanics in compari- 
son with WS model. Since the inner region of WS potential is 
nicely approximated by the power-law potential (see Fig. [TJ, 
the quantum spectra below Fermi energy for both potential 
models show fairly nice agreements in wide range of radial 
parameter a (see Fig.O. 

In Fig. [H the radial parameter a is determined so that the 
power-law potential best fit the inner region of WS potential. 
As a simple local matching, one may equate the derivatives of 
two potentials at the nuclear surface r = RA, which gives 



a = RA/2a. 



(2.3) 



Thus, the radial parameter a controls the surface diffuseness. 
For a global fitting, we take more elaborate way which min- 
imize volume integral of squared potential difference within 
nuclear radius 



— J^^ flfrr2|y„(r;A)-ywsM)| =0. (2, 

The value of a numerically obtained by Eq. ( 12.4b has an ap- 
proximately linear dependence on Ra/ci, 
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a - -0.76 + 0.69/?A/fl, 



(2.5) 



which has similar dependence on surface diffuseness a as the 
result of local fitting (I2.3l l. In Fig. |2] we use the relation 
(12.5b for the comparison of quantum spectra between WS and 
power-law potentials. Thus, we can employ this power-law 
potential model for the study of realistic nuclear shell struc- 
tures from light to heavy nuclei. 



B. Scaling properties 



FIG. 2. Single-particle spectra for spherical power-law (r") potential 
and Woods-Saxon (WS) potential as functions of radial parameter a. 
Positive and negative parity levels {n = ±) are respectively plotted 
with solid and dashed lines for power-law potential and with filled 
and open dots for WS potential. For the levels of WS model ( 12. lb . 
nuclear radius Ra is determined by Eq. ( 12.5b as a function of a. 



ing properties, which highly simplifies our semiclassical anal- 
ysis. In the following, we eliminate the constant term —W in 
Eq. ( 12.21 ) and consider the Hamiltonian for a particle of math 
m moving in the deformed power-law potential as 



u 



(2.6) 



Here, R and U are constants having dimension of length and 
energy, respectively. The dimensionless function f{6,(p) de- 
termines the shape of the equi-potential surface, and it is nor- 
malized to satisfy volume conservation condition 



(2.7) 



Under suitable scale transformation of coordinates, energy 
eigenvalue equation is transformed into a dimensionless form 



1, 



/(0,<P) 



\lf{u)=exir{u), (2.8) 



where u and e are dimensionless coordinates and energy de- 
fined by 



r E h"- 

R U mR^ 



(2.9) 



and is Laplacian in terms of the coordinate u. Since 
Eq. ( 12.81 ) does not include constants such as m, U, R and Ti, 
one can consider quantum eigenvalue problem independently 
on those values. Their absolute values are determined by fit- 
ting to the WS potential through the relation 



There is a great advantage to replace WS potential with the 
power-law potential. The power-law potential has useful scal- 
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TABLE I. Values of radial parameter a, length unit R and energy unit 
U of power-law potential l l2.6t for nuclei with mass number A. Nu- 
clear radius Ra = l.lA'/^fm, VK = 50 MeV, a = -0.76 + 1.66A'/^ 
and nucleon mass m = 938 MeV/c^ are used. 



A 


a 


R [fm] 


U [MeV] 


20 


3.75 


2.36 


3.90 


100 


6.95 


4.01 


1.33 


200 


8.95 


5.15 


0.84 



C. Semiclassical level density 

Let us consider the single-particle level density for the 
Hamiltonian ( 12. 6t . Average level density g{E) is given by 
Thomas-Fermi (TF) theory 



gTF{E) 



dpdr5{E-H{p,r)) 



2V2 /3 3\ £3 



(2.15) 



The values of a, R and U for several A are listed in TableHl 

The scaling property of the system is particularly advan- 
tageous in the analysis of classical dynamics. Since the po- 
tential is a homogeneous function of coordinates, Hamilton's 
equations of motion have the invariance under the following 
scale transformation 



{c^p,car,c«- 5fj as 



■cE. 



(2.10) 



Therefore, classical phase space structure is independent of 
energy. A phase space trajectory (ro(f),/?o(f )) energy £0 is 
transformed to a trajectory at different energy E by 



r(r) 



with 



(2.11) 



Thus we have the same set of periodic orbits in arbitrary en- 
ergy surface related through the above scale transformation. 
In the following, we set the reference energy ?XEq^U . The 
action integral along a certain periodic orbit j3 is expressed as 



pdr 



U 



SpiU)=hTp£. (2.12) 



In the last equation, we define dimensionless "scaled energy'' 
£ and "scaled period" Tp of periodic orbit j3 by 



E\ l+a 
U 



Sp{U) 



(2.13) 



The ordinary (non-scaled) period of the orbit j3 is then given 
by 



dSp{E) _ d£ 



dE 



dE 



(2.14) 



It is convenient to express periodic-orbit quantities in terms of 
£ and T in place of E and T . In HO-type potentials (a = 2), 
£ and T are proportional to ordinary energy E and period T, 
respectively. In cavities (a = 00), they are proportional to 
momentum p and orbit length L, respectively. 



which is independent of deformation under volume conserva- 
tion condition ( 12. 7t . B(s,f) represents Euler's beta function 
defined by 



B{s,t) 



\\-xy-^dx. 



By transforming energy £ to a scaled energy £, one obtains 
the scaled-energy level density 

Using ( 12.151 ), the average part is given by 

Correction to the TF density is obtained by the extended 
Thomas-Fermi (ETF) theory ill O 



1 /2m 
96^ W 
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X — / dre{E-V)- 
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For spherical case, we have the expression 



(2.18) 



(2.19) 



and the average number of levels up to energy E {£) is given 
by 

3n \ a 2J 

-^Bfl + iiV- (2.20) 
12V27r V « 2/ 

In Fig. [3] quantum-mechanically calculated coarse-grained 
level density 
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FIG. 3. Level density (upper panel) and number of levels (lower 
panel) for spherical power-law potential with a = 4.0. In the upper 
panel, solid line show the course-grained quantum level density with 
smoothing width A£ = 0.3. ETF level density given by Eq. l |2.19t 
is plotted with dashed line. In the lower panel, solid line is the sum 
of step functions l l2.22t . and dashed line is the ETF result given by 



are compared with ETF results. One sees that ETF (TF) cor- 
rectly describes the average properties of quantum results. In 
these plots, difference between ETF and TF is invisibly small. 

Next we rewrite the trace formula (11. U using scaled energy 
and scaled periods. Semiclassical formula for scaled-energy 
level density is expressed as 



(2.23) 



For an isolated orbit j3 with n repetitions, amplitude factor is 
given by the standard Gutzwiller formula SIS, 



7rr2^|det(I-M«) 



dE 
d£ 



7r,/|det(I-M"; 



(2.24) 



In the last equation, we used Eq. ( |2.14t . Mp represents the 
monodromy matrix lfl4l [Tsll . which is a linearized Poincare 
map defined by 



M« 



^K(0),/7j.(0)) 



(2.25) 



where {r±{t),p±{t)) are the coordinates and momenta per- 
pendicular to the periodic orbit j3 as functions of time f , and 
Tp is the period of the orbit. 




In a two-dimensional autonomous Hamiltonian system, 
monodromy matrix M is a (2 x 2) real and symplectic matrix. 



MJM^ =J, J 



and its eigenvalues appear in one of the following three forms: 

(a) (e",e^"): hyperbolic with no reflection 
TrM = 2 cosh M > 2 

(b) (e">-"'); elliptic 

TrM = 2 cos V, |TrM|<2 

(c) (— e",— e^"): hyperbolic with reflection 
TrM = -2 cosh M < -2 

The orbit is stable in the case (b) and otherwise unstable, and 
stability of the orbit is determined by the trace of monodromy 
matrix. The stability factor in (|2.24| i is also determined by the 
trace of monodromy matrix: 



det(I-Mp) = 2-TrM|3. 



(2.26) 



The eigenvalues of M (and therefore TrM) are independent 
of a choice of Poincare surface or a choice of canonical vari- 
ables. These eigenvalues continuously vary as deformation 
changes, and it happens that they become unity at certain val- 
ues of deformation, namely, m = in (a) or v = in (b). At 
those deformations, Poincare map acquires a new fixed point 
in the direction of eigenvector 5Zi belonging to the unit eigen- 
vector: 



M5Zi = 5Zi . 



(2.27) 



In this way, periodic orbit bifurcation occurs at TrM = 2. The 
number of new emerging orbits is dependent on the type of 
the bifurcation lfl6l . When a stable (unstable) orbit undergoes 
pitchfork bifurcation, it turns unstable (stable) and a new sta- 
ble (unstable) orbit emerges from it. When a stable orbit un- 
dergoes period-doubling bifurcation, a pair of stable and un- 
stable orbits will emerge. 

In a three-dimensional Hamiltonian system, the size of 
monodromy matrix becomes (4 x 4). Under axial symmetry, 
periodic orbits degenerate with respect to the rotation, and the 
monodromy matrix has unit eigenvalue corresponding to the 
direction of the rotation. Thus, excluding the rotational de- 
grees of freedom, stability of the orbit is described by (2 x 2) 
symmetry-reduced monodromy matrix and it has the same 
properties as in the two-dimensional case. For such degen- 
erate orbits, the trace formula is modified into what given by 
extended Gutzwiller theory [l3l. Their amplitude is pro- 
portional to the stability factor similar to that in (I2.24l i. but 
with symmetry -reduced monodromy matrix Mp . For fully de- 
generate orbits in integrable system, one can use Berry-Tabor 
formulajTll. 

In general, the stability factor | det(I — Mjg)|^'/" has strong 
dependence on the deformation parameter, and is responsi- 
ble for the sensitivity of shell structures upon deformations. 
The divergence of the Gutzwiller amplitude (12.24b based on 
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standard stationary phase method can be remedied by im- 
proved treatment of trace integral in phase space (e.g., uni- 
form approximations lfl9U2lll and improved stationary -phase 
method lllo[|22ll ) but the ampHtude may suffer strong enhance- 
ment around the bifurcation point. Since monodromy matrix 
has unit eigenvalue there, a local family of quasi-periodic or- 
bits are formed in direction of eigenvector 8Z\ belonging to 
the unit eigenvalue, and they make coherent contribution to 
the level density. Our previous studies have proved that bifur- 
cations of short periodic orbits play significant role in emer- 
gence of deformed shell structures 1191 [lOl I22I - I24I1 . 

One should however note that the above enhancement is not 
always found for every bifurcations. The significance of bifur- 
cation depends on the normal form parameters which describe 
nonlinear dynamics around the periodic orbit at the bifurca- 
tion points. In Ref. ll25ll . uniform approximation remedies 
the divergence problems which one encounters at bifurcation 
points in standard stationary phase method, but the obtained 
amplitude show no enhancement around there. In Ref. ll24ll . 
we found very strong enhancement of amplitude around the 
bifurcation point for one certain orbit, but the same type of 
bifurcation in another orbit shows no enhancement. In our 
experience, significant growth of shell effect at a certain de- 
formation is related with bifurcations of simple short periodic 
orbits. 



D. Fourier transformation technique 

Fourier transformation technique is especially useful in 
studying classical-quantum correspondence in the system 
with scale invariance. Let us consider the Fourier transform 
of scaled-energy level density 

F{t) = j g{£)e''^e-^-^f^^~d8. (2.28) 

In the integrand, Gaussian is included in order to exclude 
the level density at high energy S> 1 where single-particle 
spectra are not obtained numerically in a good precision. 

By inserting the quantum level density g{£) = 8{£ — £„) 
into Eq. ( 12.28b . one obtains 

which can be easily calculated using quantum mechanically 
calculated energy eigenvalues {£„]. On the other hand, by 
inserting the semiclassical level density ( 12.23b . one formally 
has the expression 

F^i(T) =F(T) + 7r£e"^"/^/2A„/5(-'^.)5r(T-«T^). (2.30) 
Here, 5y{z) represents a normalized Gaussian with width 7 

1 -4 

5y{z) = ^^e ^, (2.31) 

which coincides with Dirac's delta function in the limit 7^-0. 
Thus, F{t) is a function possessing successive peaks at the 
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FIG. 4. Shell energies calculated for power-law potential models 
with a = 4.0, 5.0 and 8.0, plotted as functions of cubic root of par- 
ticle number A', which is either neutron or proton number, taking 
account of the spin degeneracy factor. 

scaled periods of periodic orbits T = iiTp. Therefore, by cal- 
culating Fourier transform of quantum mechanical level den- 
sity, one can extract information on the significance of each 
periodic orbits contributing to semiclassical level density. 

The parameter 7 implies the resolution of the periodic or- 
bit in the Fourier transform. For a better resolution, larger 
number of good quantum energy levels (up to £max 2/7) are 
required. 



III. SPHERICAL POWER-LAW POTENTIALS 
A. Quantum spectra and shell structures 

Let us first examine the quantum spectra in spherical power- 
law potentials. In case of a = 2 (isotropic harmonic oscilla- 
tor), energy levels with different angular momentum degen- 
erate due to dynamical SU(3) symmetry. With increasing 
a, these degeneracies are resolved and shell effect becomes 
weaker. However, new types of prominent shell structures will 
emerge at certain values of a. The shell structure, a fluctua- 
tion in single-particle spectrum, brings about a fluctuation in 
energy of nuclei as functions of constituent nucleon numbers. 
This fluctuation part, which we call shell energy, is calculated 
by removing the smooth part from a sum of single-particle 
energies by means of Strutinsky method 1261 127I1 . 

Figure |4] displays shell energies as functions of particle 
(neutron or proton) number for several values of a. For 
a = 4.0, one finds subshell structures with 3 bumps per each 
major shell which are more evident in large regions. For 
a = 5.0, the subshell structure changes into that with 2 bumps 
per each major shefl. For a = 8.0, the magnitude of major 
shell effect is considerably enhanced. The last one may in- 
dicate an enhancement of shell effect around medium-mass 
region. In the following subsections, we will show that the 
above unique shell structures can be explained as the effect of 
different types of periodic orbit bifurcations which take places 
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FIG. 5. Some short periodic orbits {k,n) in spherical power-law po- 
tentials. They emerge from the circular orbit (shown in broken lines) 
via period «-upling bifurcations. 
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around the above values of a. 



B. Classical periodic orbits 

In spherical power-law potential model, several simple an- 
alytic descriptions for the properties of the periodic orbits are 
available. Taking the orbits in {x,y) plane and putting the z- 
component of angular momentum to /- = K, 2-dimensional 
effective Hamiltonian in polar coordinates is written as 



H=-pj + V,s{r,K), V,ff{r,K)=r" + -^ 



(3.1) 



Circular orbit r(f ) = (denoted by C) satisfies the condition 



0, 



from which one obtains, for energy E, 

l/a 



2E 
l + a 



and the angular frequency 



(Oc 



K 

72 



2E 
l + a 



I 1 



(3.2) 



(3.3) 



(3.4) 



The circular orbit is stable, and r(f ) of the orbits in vicinities 
of the circular orbit oscillate around with angular frequency 



eff 



2E 



a 



(3.5) 



Bifurcations occur when the ratio of those two frequencies cOc 
and Qc becomes rational, namely. 
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FIG. 6. Modulus of Fourier transform of quantum level density 
l |2.29| l plotted as functions of scaled period t for a = 4.0, 5.0 and 
8.0. 



for period m-upling bifurcation. Here, a new orbit which os- 
cillate n times in radial direction when it rotates m times along 
the orbit C emerges from m-th repetition of orbit C. The values 
of a at such bifurcations are given by 



a = 



(3.7) 
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Figure |5] shows some periodic orbits («,m) emerge from the 
circular orbit via period multipling bifurcations of circular or- 
bit C. 

In Fig.|6]we show Fourier amplitudes of quantum level den- 
sity ( |2.29t plotted as functions of scaled period t for several 
values of a. In the upper panel (a = 4.0), one sees a large 
peak at T ~ 17. This peak corresponds to 3rd repetitions of 
circular orbit, 3C, as well as orbit (7,3) bifurcated from 3C at 
a = 31/9 = 3.44. Since these orbits have periods about three 
times as long as primitive circular orbit, the interference of 
their contribution to the semiclassical level density will show a 
shell structure such that a major shell is modulated by a small 
oscillations having triple frequency. Thus, the appearance of 
subshell structure found in the upper panel of Fig. 2] is un- 
derstood as a result of bifurcation enhancement effect of orbit 
(7,3). In the middle panel, the highest peak at T = 12 corre- 
sponds to the 2nd repetitions of circular orbit, 2C, as well as 
orbit (5,2) bifurcated from 2C at a = 4. In the lower panel, 
peak at T ^ 7 corresponds to the primitive circular orbit C 
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FIG. 7. Gray-scale plot of Fourier transform of quantum level density 
l |2.29l l as function of radial parameter a and scaled period T. Mod- 
ulus of the Fourier transform |F(t)| has large value at dark region. 
Scaled periods of classical periodic orbits (a) are also shown with 
lines as functions of a. Bifurcation points are indicated by open cir- 
cles. 



FIG. 8. Poincare surface of section {{x,px)\z = 0} for meridian- 
plane orbits in spheroidal-shape potentials with axis ratio rj = V3 
(8 ~ 0.49) and with several values of radial parameter a. The ori- 
gin (jc = px = 0) corresponds to the orbit Z, and the outer boundary 
corresponds to the orbit X. 



and orbit (3,1) bifurcated from C at a = 7. Shell structures 
found in the middle and lower panels of Fig. |4] can be under- 
stood, just in the same way as upper panel, as the bifurcation 
enhancement effect of orbit (5,2) and (3,1), respectively. 

Figure |2] shows a two-dimensional gray-scale plot of quan- 
tum mechanical Fourier amplitude |F(T;a)| as a function of 
radial parameter a and scaled period T. Scaled periods of 
classical periodic orbits are also drawn as functions of a. We 
see nice agreements between quantum Fourier transforms and 
classical periodic orbits. The Fourier peaks associated with 
periodic orbits are strongly enhanced around their bifurcation 
points indicated by open circles. This clearly illustrates the 
significance of periodic-orbit bifurcations to the enhancement 
of shell effect. 



IV. SPHEROIDAL DEFORMATIONS 

A. Shape parametrization and quantum spectra 

Axially-symmetric anisotropic harmonic oscillator poten- 
tial system is integrable, and it has spheroidal equipotential 
surface. It is known that spheroidal deformed cavity (infinite 
well potential) system is also integrable. For spheroidal de- 



formation, shape function is written as 

/(0) = 



sin^e 



cos^ 9 



-1-1/2 



[R^/RqY {RjRoy 



(4.1) 



where R- and R± represent lengths of semiaxes of the spheroid 
which are parallel and perpendicular to the symmetry axis (z- 
axis), respectively. Taking account of the volume conservation 
condition R\R^ = Rq, we define the deformation parameter 5 
as 



Ri 



Roe 



-5/3 



R: = Roe 



25/3 



(4.2) 



It is related with the axis ratio Tj = Rz/Rj- by Tj = e^. Spher- 
ical shape 77 = 1 corresponds to 5 = and prolate/oblate su- 
perdeformed shapes 77 = 2^' corresponds to 5 = ±log2 w 
±0.69. The system with spheroidal power-low potential is 
non-integrable except for two limits a = 2 (HO) and a = 00 
(cavity). In Fig. [8] we show Poincare surface of section for 
a = 2,5,20 and 00, each with t] = ^3 (5 « 0.49). With 
increasing a > 2, some complex structures emerge in the 
Poincare plots and it becomes most chaotic around a ^ 20, 
then turns into simpler structure for extremely large a. 
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FIG. 9. Single-particle level diagram for Hamiltonian ( 12. 6b with 
spheroidal deformation J4. It . Scaled-energy levels are plotted as 
functions of deformation parameter 8 defined by Eq. i4.2i . Solid 
and broken lines represents levels with even and odd parities, respec- 
tively. 



Figure |9] shows the single-particle spectra as functions of 
spheroidal deformation parameter 5. The value of radial pa- 
rameter is put a = 5.0, corresponding to medium-mass nu- 
clei. The degeneracies of levels at spherical shape are resolved 
and shell structure changes as varying deformation. The level 
diagram is similar to what obtained for MO or WS models 
without spin-orbit coupling. One of its characteristic features 
in comparison to HO model is the asymmetry of deformed 
shell structures in prolate and oblate sides. This asymmetry 
becomes more pronounced for larger a, and it might be re- 
garded as the origin of prolate-shape dominance in nuclear 
ground-state deformations. We shall discuss the semiclassical 



O=5,ft=0.3 



a=5,ft=a4 




■i 



FIG. 11. Same as Fig.|9]but as functions of quadrupole deformation 
parameter j32. 



origin of the above asymmetry in the following subsections. 

In order to see the dependence on shape parametriza- 
tion, we also calculated the deformed quantum spectra for 
quadrupole deformation, which might be more popular in ear- 
lier studies, 



m 



(4.3) 



The factor in the denominator arranges the conservation of 
volume surrounded by equipotential surface. Figure [TO] show 
Poincare surface of section for quadrupole deformations j82 = 
0.3 and 0.4 with a = 5.0. Comparing with right upper panel 
in Fig. [8] one will see that the particle motions in quadrupole 
potential are more chaotic than those in spheroidal potential. 

Figure[TT]shows the level diagram for quadrupole deforma- 
tion. Although the properties of classical motions are quite 
different from those in spheroidal potential, the deformed 
shell structures are very similar to each other. Thus, the above 
difference of shape parametrization does not cause a serious 
difference in the gross shell structures at normal deformations. 
Notable effects of chaoticity in quadrupole potential can only 
be seen in the strong level repulsions at large deformations 

^ 0.3. Therefore, we shall only consider the spheroidal 
deformation in the following analysis. 



B. Prolate-oblate asymmetry in deformed shell structures 



FIG. 10. Same as Fig.[8]but for quadrupole deformations P2 = 0-3 
and 0.4 with a = 5.0. 



Let us examine the properties of deformed shell structures 
in normal deformation region ( 1 5 1 < 0.3). As shown in Figs. |9] 
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and compare the energies in prolate and oblate sides at each 
local minima. Here, we assume the same single-particle spec- 
tra for neutrons and protons and only consider N = Z even- 
even nuclei for simplicity. The sum of single-particle energies 
for nucleus of mass number A is given by 

n 

£.vp(A) =4^e/, A=N + Z^4n. (4.5) 

i=l 

Using the Strutinsky method, the above energy can be decom- 
posed into smooth part Esp{A) and oscillating part 5E{A). As 
in usual, we can expect that the above oscillating part repre- 
sents the correct quantum shell effect of many-body system. 
In Strutinsky's shell correction method, the smooth part is re- 
placed with the phenomenological liquid drop model (LDM) 
energy to get the full many-body energy, but here we try to 
extract the smooth part also from the single-particle energies. 
In mean-field approximation, single particle Hamiltonian is 
written as 



ii = i+r 



(4.6) 



FIG. 12. Fanning of single-particle levels Ih (dashed lines) and li 
(solid lines) for a = 5.0 (lower panel) and a = 1 . 1 (upper panel). 



and[TT] single-particle spectra in a potential with sharp surface 
show prolate-oblate asymmetry. Hamamoto and Mottelsonl3i| 
paid attention on the different ways of level spreading (fan- 
ning) in oblate and prolate sides; level spreading is consid- 
erably suppressed in the oblate side as compared to the pro- 
late side. Due to that suppression of level spreading, shell 
structures in the oblate shapes are similar to that of spher- 
ical shape and system has smaller chance to gain shell en- 
ergy by means of oblate deformation. This may explain the 
feature of prolate-shape dominance. They have shown that 
the above asymmetric way of level fanning can be understood 
from the different roles of the interaction between low-K lev- 
els for level repulsions in oblate and prolate sides in a potential 
with sharp surface. It clearly explains the fact that the asym- 
metry becomes morepronounced for heavier nuclei e.g., in 
Woods-Saxon model@]. The same kind of asymmetry is also 
found in the spectrum of Nilsson model. 

In spheroidal power-law potential model, the asymmetry 
in level fanning becomes more pronounced for larger a as 
expected. In the lower panel of Fig. [121 spreading of some 
nl levels (« and / represent principal and azimuthal quantum 
numbers, respectively) is illustrated for a = 5.0. One sees that 
the level spreading is considerably suppressed in oblate side 
as in cavity potential. 

It is interesting to note that if we take the radial parameter 
a <2 (although it does not correspond to actual nuclear situ- 
ations), the way of level fanning becomes just opposite to the 
case of a > 2. As one see in the upper panel of Fig. [T2] l. level 
spreading is suppressed in the prolate side. One may ask if 
that should cause oblate-shape dominance. 

Following the analysis in Ref. IH, we calculate the defor- 
mation energy 



£def(A,5)=£(A,5)-£(A,0) 



where f and F represent kinetic energy and mean-field poten- 
tial, respectively, and F is currently given by the power-law 
potential. In this case, by the use of Virial theorem, average 
of t and F are in ratio 2(f) = a(r), and one obtains 



a 



it = 



a + 2 



{h] 



a + 2 



(4.7) 



Therefore, the smooth (average) part of the A-body energy is 
given approximately by 



E{A) 




a 



a 



(4.8) 



This expression will be valid for many-body systems interact- 
ing with 2-body interaction. Thus, we calculate the A-body 
energy by 



£(A) 



a + 1 



a 



E,p{A) + 5E{A) 



(4.9) 



(4.4) 



Figure [13] compares the local minima of deformation en- 
ergies ( 14.41 ) in prolate and oblate sides. At the HO value 
(a — 2.0), prolate and oblate deformed shell structures are 
symmetric and the deformation energies are comparable with 
each other. For a > 2 (two bottom panels), the deformation 
energies in prolate side become considerably lower than in 
oblate side as the radial parameter a becomes larger The 
power-law potential model thus reproduce correctly the fea- 
ture of prolate-shape dominance in nuclear deformation. 

For a < 2, as shown in the top panel of Fig. [T3] we find 
no indication of oblate-shape dominance in spite of the fea- 
ture of level fanning shown in the upper panel of Fig. [12] 
Some nuclei have lowest energies at oblate shapes 5 ^ —0.3, 
but the difference in energies between prolate and oblate min- 
ima are generally small. Therefore, one cannot fully explain 
the prolate(oblate)-shape dominance only by the ways of level 
fanning. 
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FIG. 13. Deformation energies i4Al of prolate and oblate states at 
each local minima. Deformation energy £a(5) is calculated as func- 
tion of deformation 5 and its minimum values for prolate and oblate 
sides are plotted with filled and open circles, respectively. 



In order to analyze shape stability, we define shell- 
deformation energy using the smooth part of energy at spher- 
ical shape as reference. 



AE{A,5)=E{A,5)-E{A,Q) 



(4.10) 



with ( 1481 ) and i43[ . (Note that the 2nd term of ( 14.10b is not 
£■(^,5), so that AE contains smooth part of deformation en- 
ergy.) Figure [T4l shows contour plots of AE for a = 1.1 and 
5.0 as functions of deformation 8 and mass number A. They 
show some deep minima on 5 = line at values of A cor- 
responding to spherical magic numbers. The valley lines run 
through these minima and the deformation energy minima dis- 
tribute along these valleys. For a = 5.0, the valley lines have 
large slopes in the prolate side and deep energy minima are 
formed at 5 ^ 0.2 between spherical magics, while the valley 
lines are almost flat in the oblate side. For a = 1.1, the valley 
lines have larger slope in the oblate side, but the slope in pro- 
late side is not as small as in oblate side for a = 5.0 and the 
deformation energy minima distribute mainly along the val- 
ley lines in the prolate side. One can find rather deep energy 
minima at 5 '--^ —0.3 for particle numbers between spherical 
magics, but the energy difference between oblate and prolate 
local minima are generally small. Thus, for an understanding 
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FIG. 14. Contour plot of the A-body shell-deformation energy loot 
in deformation-mass number plane {5,A^^^) for a = l.I (upper 
panel) and a = 5.0 (lower panel). Solid and dashed contour lines 
represent negative and positive values, respectively. Dots represent 
values of deformation parameter at absolute energy minima for each 
A. 



of prolate-shape dominance, it is critical to explain the asym- 
metric behavior of the slopes of the energy valleys. 

Since shell energy takes deep negative value when the 
single-particle level density at the Fermi energy is low, let us 
investigate the coarse-grained single-particle level density as 
functions of energy and deformation. Figure [15] shows the 
oscillating part of course-grained single-particle level density 
for a = 1 . 1 and 5 .0 plotted as functions of deformation 8 and 
scaled energy £. They show regular ridge-valley structures 
similar to shell energy. Therefore, for an understanding of 
prolate-shape dominance, it is essential to investigate the ori- 
gin of the above ridge-valley structures in deformed single- 
particle level density. 
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FIG. 15. Contour plot of oscillating level density for radial parameter 
a = 1.1 (upper panel) and 5.0 (lower panel) in deformation-energy 
plane {S,£). For calculation of coarse-grained level density, smooth- 
ing width A£ = 3.0 is taken. Solid and dashed contour lines represent 
negative and positive values, respectively. Thick solid lines represent 
constant-action lines l |4.13t for the bridge orbits M(l,l). 

In the following subsections, we will show that the above 
ridge-valley structure can be explained in connection with 
classical periodic orbits using semiclassical periodic-orbit the- 
ory. The thick solid curves in Fig. [15] represent the semiclas- 
sical prediction of the valley lines which will be discussed in 
Sect. UVDl 



C. Periodic orbits in spheroidal potential 

In order to examine the semiclassical origin of the above 
asymmetry in deformed shell structure using periodic orbit 
theory, we first consider the properties of classical periodic 



orbits in spheroidal power-law potential and their bifurca- 
tions. For spherical potential, all the periodic orbits are pla- 
nar and degenerate with respect to rotations. The degener- 
acy parameter is /C = 3 for generic orbits, and K, = 2 for di- 
ameter and circular orbits which are transformed onto them- 
selves by one of the rotations. If the spheroidal deformation 
is added to the potential, generic planar orbits bifurcate into 
two branches: One is the orbit in equatorial plane and the 
other is the orbit in meridian plane (plane containing sym- 
metry axis). All but two exceptional orbits degenerate with 
respect to the rotation about symmetry axis, and the degen- 
eracy parameter is /C = 1 . The diametric orbit bifurcate into 
degenerate equatorial diameter (/C = 1) and isolated diametric 
orbit along symmetry-axis (JC ~ 0). The circular orbit bifur- 
cate into isolated equatorial circular orbit (/C = 0) and oval 
orbit in meridian plane {JC = I). With increasing deforma- 
tion towards prolate side {8 > 0), the equatorial orbits un- 
dergo successive period-multipling bifurcations and new 3D 
orbits emerge. In the oblate side, diametric orbit along the 
symmetry-axis undergoes successive period-multipling bifur- 
cations and generate new meridian-plane orbits. These new- 
born 3D and meridian-plane orbits have hyperbolic caustics 
and are sometimes called "hyperbolic orbits". 

It is very interesting to note that the above new-born hyper- 
bolic orbits from equatorial orbits are distorted towards sym- 
metry axis as increasing deformation and finally submerge 
into diametric orbit along symmetry-axis. (Some 3D orbits 
submerge into other hyperbolic orbits before submerging into 
symmetry-axis orbit at last.) In this way, the hyperbolic or- 
bits make bridges between the equatorial and symmetry-axis 
orbits, and we shall call those hyperbolic orbits as "bridge 
orbits" ll28ll29ll . As increasing 8, periods of equatorial orbit 
decrease while that of symmetry-axis orbit increases. At each 
crossing point of the periods (or actions) of repeated equato- 
rial and symmetry-axis orbits, bridge orbit exist to intervene 
them. 

Accordingly, we shall classify periodic orbits in spheroidal 
power-law potential into the following 4 groups: 

i) Isolated orbits {JC = 0): This group consists of the dia- 
metric orbit along symmetry axis (z-axis), denoted "Z", 
and the circular orbit in equatorial plane, denoted "EC". 
Orbit EC is stable both in the prolate and oblate sides, 
whose repeated period-multipling bifurcations gener- 
ate 3D bridge orbits. Orbit Z is stable in the oblate 
side and undergoes successive period multipling bifur- 
cations, while its stability alternate in the prolate side 
with repeated bifurcations which absorb bridge orbits. 

ii) Equatorial-plane orbits {JC — 1): This corresponds to 
the equatorial-plane branch of the deformation-induced 
bifurcation. They have the same shapes as those in 
spherical potential shown in Fig. |5] They are denoted 
E(^,m), where k is the number of vertices (corners) and 
m is the number of rotation. Diametric orbit is espe- 
cially denoted "X" (which include the orbits along x 
axis). 

iii) Meridian-plane orbits {JC = 1): This corresponds to the 
meridian-plane branch of the deformation induced bi- 
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FIG. 16. Bifurcation diagram for M(l,l) bridge orbit between X 
and Z orbits for a = 3.0. Traces of symmetry-reduced monodromy 
matrix M are plotted as functions of deformation parameter d. Bi- 
furcation points (TrM = 2) are indicated by open circles. Shapes of 
the periodic orbits as well as equipotential surface are also shown. 

furcation. They survive up to any large deformation, 
keeping their original geometries. 

iv) Bridge orbits {K. — 1); These orbits emerge from the 
bifurcations of equatorial orbits. Meridian-plane orbits 
emerge from diametric orbits and submerge into rep- 
etitions of orbit Z. Non-planar 3D orbits emerge from 
non-diametric equatorial orbits, and they also submerge 
into the orbit Z. Some of them submerge into other 
bridge orbit before submerging into Z. The meridian- 
plane bridge intervening mX (m-th repetition of X) and 
nL («-th repetition of Z) is denoted M(m,n). Except 
for M(l,l) bridge, a pair of stable and unstable bridge 
orbits emerge, and are denoted M(ot,«)j and M(m,n)„, 
respectively. 3D bridge B(m,m,n) emerge via pitchfork 
bifurcation of equatorial mEC (m-th repetitions of EC), 
and submerge into M(m,«) orbit before finally sub- 
merge into «Z. The other 3D bridges intervening equa- 
torial E(A:,m) and «Z emerges as a stable and unstable 
pair, and denoted as B(A:,;7i, «)(,„. As increasing defor- 
mation, they first submerge into 3D bridge B(OT,OT,n), 
which will submerge into M(OT,n) and finally into nZ. 

Figure [16] shows bifurcation diagram for bridge orbit 
M( 1,1) for a = 3.0. The traces of (2 x 2) symmetry -reduced 
monodromy matrices for relevant periodic orbits are plotted 
as functions of deformation parameter 5. The JC — \ family 
of equatorial diameter orbit X undergoes pitchfork bifurcation 
at 5 = —0.34 and a family of oval-shape meridian-plane orbit 
]V1(1,1) emerge. In the limit 5 0, the shape of M(l,l) orbit 
approaches circle, and it associates with equatorial circular 
orbit EC to form JC = 2 family. At 5 > it bifurcates into 
equatorial EC and meridian M(l,l) family again. The merid- 
ian branch submerge into the orbit Z at 5 = 0.34 via pitchfork 
bifurcation. Thus, the orbits M(l,l) make bridge between two 
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FIG. 17. Modulus of Fourier transform of the level density is shown 
by gray-scale plot as a function of deformation parameter S and 
scaled period T. Scaled periods of classical periodic orbits are also 
displayed. 



diametric orbits X and Z. 

In the HO limit, a ^ 2, the bridge shrinks to a crossing 
point of X and Z orbits and can exist only at 5 = (spherical 
shape), where they altogether form a degenerate /C = 2 family. 
With increasing a, the bridge orbit exist in wider range of de- 
formation over the crossing point. In the cavity limit, a oo, 
this orbit approaches so-called creeping orbit or whispering 
gallery orbit, which runs along the boundary. 



D. Semiclassical origin of prolate-oblate asymmetry 

To see the effect of the above bifurcation to the shell struc- 
ture, we calculate Fourier transform of level density (|2.29t 
with the obtained quantum spectra. In Fig. [TT] modulus of 
Fourier transform |F(t,5)| is shown in gray-scale plot as a 
function of deformation 8 and scaled period T. Scaled peri- 
ods of classical periodic orbits Zp (5) are also drawn by lines. 
One sees nice correspondence between quantum Fourier am- 
plitude and classical periodic orbits. Particularly, one can find 



significant peaks along the bridge orbit M( 1,1), which indicate 
that the shell structure in normal deformation region is mainly 
determined by the contribution of this bridge orbit. 

Let us assume that a contribution of single orbit (or degen- 
erate family) j3 dominate the periodic orbit sum, namely, 

5g{£)^Apcosi£Tp-^Vp). (4.11) 

Then, the valley lines of level density should run along the 
curves where above cosine function takes the minimum value 
— 1 , namely, 

£Tp~^Vp^{2n + l)K, «==0,1,2,---. (4.12) 
In Fig. [T3] we plot this constant-action curves 

£ = - T^r^ (4-13) 

for bridge orbit M(l,l). We see that the constant-action lines 
of the bridge orbit nicely explain the ridge-valley structure in 
quantum level density. Properties of deformed shell energy 
minima in Fig. [14] are thus understood as the effect of bridge 
orbit contribution. For a > 2, bridge orbits appear upward 
from the crossing point of two diametric orbits X and Z in 
(5, t) plane. Due to the volume conservation condition, orbit 
Z has larger slope than X orbits. Therefore, the bridge be- 
tween X and Z orbits have larger slope in prolate side than in 
oblate side. This explains the origin of prolate-oblate asym- 
metry in deformed shell structure. For a < 2, bridge orbit 
appear in opposite side of the crossing point and its slope be- 
comes larger in oblate side. This also explain the nature of 
valley lines in level density and shell energy for a = 1.1 in 
Figs.[T5]and[ni 



E. Superdeformed shell structures 

In harmonic oscillator potential model, one sees simultane- 
ous degeneracy of many energy levels at rational axis ratios. 
Superdeformation is formed under the strong level bunching 
at axis ratio 2:1. Search of much larger deformation originated 
from the shell structure at axis ratio 3:1 (sometimes referred to 
as hyperdeformation) has also been a challenging experimen- 
tal and theoretical problem. In spheroidal cavity system, one 
also sees superdeformed shell structure and it is intimately re- 
lated with emergence of 3D orbits which turns twice around 
symmetry axis while it oscillate once along symmetry axis, 
just as the degenerate 3D orbits in 2:1 axially-deformed HO 
potential^. One may expect to have a general semiclassical 
understandings for the origin of superdeformed shell struc- 
tures in realistic nuclear potentials by connecting the above 
two limits with power-law potential model. 

Figure [18] shows the oscillating part of the coarse-grained 
level density for radial parameter a = 5.0 and deformation 5 
around superdeformed region. It clearly show that new regu- 
larities in shell structure are formed at superdeformed region. 
The valley lines are up-going till 5 ^ 0.5, and they bent down 
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FIG. 18. Contour plot of oscillating level density for radial parame- 
ter a = 5.0 around prolate superdeformed region. Smoothing width 
A£ = 0.2 is taken. Solid and dashed contour lines represent negative 
and positive values, respectively. Thick dotted, dashed and solid lines 
represent constant action lines l l4.13t of periodic orbits 2X, M(2,l) 
and B(2,2,l), respectively. 

around 5 ^ 0.6. One sees another deep minima at 5 > 0.7. 
Let us examine their semiclassical origins. 

For a > 2, one finds bridge orbits M(2,l) which intervenes 
orbits 2X (2nd repetitions of X) and Z. Figure [19] is the bi- 
furcation diagram for the orbits relevant to this bifurcation, 
calculated for a = 3.0. The orbit X undergoes period dou- 
bling bifurcation at 5 = 0.55 and emerge a pair of bridge 
orbits M(2,l)i (stable) and M(2, 1)„ (unstable). They have 
shapes of boomerang and butterfly as shown in Fig. [19] As in- 
creasing 5, those orbits are distorted towards z axis and finally 
submerge into the orbit Z at different values of 5 via pitchfork 
bifurcations. 

For larger a, various equatorial orbits appear as shown in 
Fig. [5] and they also undergo bifurcations by imposing de- 
formation. Each of those bifurcation will generate a pair of 
3D bridge orbits, which are also distorted towards symmetry 
axis by increasing 5 and finally submerge into Z. Figure [20] 
shows some 3D bridge orbits important for superdeformed 
shell sti'uctures for a = 5.0. Equatorial circular orbit EC un- 
dergoes period-doubling bifurcation which is peculiar to 3D 
system and generate 3D bridge orbit B(2,2,l). Equatorial or- 
bit E(5,2) undergoes non-generic period doubling bifurcation 
and a pair of 3D bride orbits (5,2, l)^,,, emerge. All the above 
3D orbits finally submerge into Z orbit by increasing defor- 
mation 5. See Appendix [A] for detailed description of these 
3D bridge orbit bifurcations. 

Figure |2T] shows Fourier transform of scaled-energy level 
density for a = 5.0 around superdeformed region. The scaled 
period of classical periodic orbits are also drawn with lines. 
The Fourier amplitude shows remarkable enhancement along 
the bridge orbits M(2,l) and B(5,2,l), indicating their signifi- 
cant roles in superdeformed shell structures. 
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0.5 0.6 0.7 0.8 



FIG. 19. Same as Fig. [16] but for M(2,l) bridge orbit between 2X 
and Z, "2X" represents 2nd repetition of X. The radial parameter 
a = 3.0 is used. A pair of bridge orbits emerge at 5 = 0.55 via 
period-doubling bifurcation. The unstable branch M(2, l)i, and sta- 
ble branch M(2, l),s submerge into Z at 5 = 0.80 0.97, respectively, 
via pitchfork bifurcations. 





(5,2,l)ii 



(2,2,1) 

FIG. 20. 3D orbits responsible for superdeformed shell structure at 
S ~ 0.7 and a = 5. 3D plots and projections on {x,y), (x,z), (v,z) 
planes are shown as well as equipotential surfaces. B(5,2, l)j and 
B(5,2, l)i, are pair of stable and unstable 3D orbits emerge from 
equatorial orbit E(5,2). B(2,2, 1) emerge from 2nd repetition of equa- 
torial circular orbit, 2EC. 
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a=5.0 . 





0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
5 

FIG. 21. Gray-scale plot of Fourier transform of the quantum level 
density j2.29t for radial parameter a = 5.0 as function of deforma- 
tion 5 and scaled period T. The modulus of the Fourier transform has 
large value at dark region. The scaled periods of classical periodic 
orbits are displayed with lines. Their bifurcation points are indicated 
by open circles, left panel represent the moduli of Fourier transform 
|F(T)j as function of T and deformation 5. The scaled period of 
classical periodic orbits responsible for superdeformed shell struc- 
ture is plotted in the right panel. 



The constant action lines (|4TT3] | for M(2,l) and B(2,2,l) 
are shown in Fig. [18] with thick solid and broken lines. They 
perfectly explain the ridge-valley structures of quantum level 
densities. It shows the significant roles of bifurcations of 
M(2,l) and («,2, 1) orbits for enhanced shell effect at 5 ^ 0.5 
and 0.7, respectively. 

M(2,l) and (2,2,1) orbits shrink to the crossing point of 
2X(2EC) and Z orbits in the HO limit, a — > 2, and turns a /C = 
4 degenerate family. As increasing a, the deformation range 
in which bridge orbit can exist becomes wider. Therefore, the 
bifurcation deformation of the orbits 2X and 2EC becomes 
smaller as increasing a and the effect of these orbits take place 
at smaller deformation. This may explain the experimental 
fact that deformation of superdeformed band is smaller for 
heavier nuclei; e.g., ^ 0.6 for Dy region and ^2 = 0.4 ~ 0.5 



forHgregion|l3a[31| 



V. SUMMARY 

We have made a semiclassical analysis of radial power-law 
potential with spheroidal deformation. We have shown that 
bridge orbits mediating equatorial and symmetry-axis orbits 
play significant role in normal and superdeformed shell struc- 
tures. Particularly, prolate-oblate asymmetry of deformed 
shell structures, which is responsible for the prolate domi- 
nance in nuclear deformations, is clearly understood as the 
asymmetric slopes of bridge orbits in {8, t) plane. This asym- 
metry grows as increasing radial parameter a, namely, as 
larger mass number A, which explain the fact that the prolate 
dominance is more remarkable in heavier nuclei. 

Establishing semiclassical method to evaluate contribution 
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of classical periodic orbit bifurcation is also an important sub- 
ject. Some preliminary results for spherical power-law po- 
tential using improved stationary phase method have been re- 
ported in II32I1 . 

Since nuclear mean field has strong spin-orbit coupling, it 
is crucial to take account of its effect to analyze realistic nu- 
clei. This will be done in Part. II, where we will expand the 
model Hamiltonian to incorporate spin-orbit coupling and dis- 
cuss the nuclear problems which are closely related with spin 
degree of freedom. 
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eigenvalue generally has off-diagonal element v: 











1 








M 



For finite v, there is only one eigenvector belonging to the 
unit eigenvalue, corresponding to the direction of symmet- 
ric rotation. This off-diagonal element varies as a function 
of deformation, and vanishes at the bifurcation point, where 
M acquires a new eigenvector perpendicular to the former 
one. Here, the symmetry-reduced monodromy matrix M gen- 
erally does not have unit eigenvalues. This is what occur in 
the case of a 3D orbit bifurcations in axially symmetric poten- 
tial, which is not detected from the trace of symmetry-reduced 
monodromy matrix. 



Appendix A: Bifurcations of 3D bridge orbits 

For a periodic orbit in 2D autonomous Hamiltonian sys- 
tem, one can examine its bifurcation scenario by evaluating 
the trace of (2 x 2) monodromy matrix as a function of con- 
trol parameter such as deformation, strength of external field, 
or energy. 3D orbits in axially symmetric potential has (2x2) 
symmetry -reduced monodromy matrix, but the ignored degree 
of freedom corresponding to symmetric rotation also play a 
role in bifurcation. 

Figure |22] shows bifurcation diagram of periodic orbits re- 
sponsible for superdeformed shell structures for a = 5.0. The 
orbit X undergoes period-doubling bifurcation at 5 = 0.42 
and generate a pair of bridges M(2, 1)„ and M(2, 1),, which 
submerge into the orbit Z at 5 = 0.78 and 0.94, respectively. 
Equatorial circular orbit EC undergoes period-doubling bifur- 
cation at 5 = 0.6 and generate 3D bridge B(2,2,l), which sub- 
merge into M(2, l)s at 5 = 0.77 before finally submerge into 
Z. Since EC is isolated, monodromy matrix has size (4 x 4) 
and its four eigenvalues consists of two conjugate/reciprocal 
pairs. One of them are e^" ' , which represent stability against 
displacement in equatorial plane, whose values are indepen- 
dent of deformation 5 (2EC(1) in Fig. |22] |. The other pair 
e±»'j^ which represent stability against displacement towards 
off-planar direction, change their values as function of de- 
formation (2EC(2) in lower panel of Fig. |22] |. Bifurcation 
occurs when the latter eigenvalues become unity (vj = 0). 
The monodromy matrix of bridge B(2,2,l) has eigenvalues 
(e"'^,e~"'^, 1, 1) at its birth, and the first two eigenvalues 
change as increasing deformation. Therefore, the bifurcation 
point does not correspond to TrM = 2 for this bifurcation. 
The orbit B(2,2,l) submerge into M(2,l) at 5 = 0.77. This bi- 
furcation point does not correspond to TrM = 2, either Here, 
with decreasing 5, the mother orbit M(2, 1 )i pushes out a new 
orbit B(2,2,l) in direction of the eigenvector of M belonging 
to one of the unit eigenvalues (other than which corresponding 
to the rotation about symmetry axis). 

In general, real symplectic matrix M can be transformed 
into Jordan canonical form by a suitable orthogonal trans- 
formation, and its (2 x 2) sub-block associated with the unit 
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FIG. 22. Bifurcation diagram for periodic orbits around superde- 
formed region for a = 5.0. Values of traces of the symmetry-reduced 
monodromy matrices are plotted as functions of deformation param- 
eter 5. The lower panel is the magnified plot of the rectangular region 
indicated in the upper panel. For the orbit 2EC, traces of two (2 x 2) 
sub-blocks in total (4 x 4) monodromy matrix, denoted by 2EC(1) 
and 2EC(2), are plotted (see text). In the lower panel, TrM for orbits 
marked '*' and '**' are plotted in different scales indicated on the 
right vertical axis. 
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